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Abstract 

In [1] we defined a function (we called S) associated to a rational second order 
ordinary differential equation (rational 20DE) that is linked to the search of an 



^Nq . integrating factor. In this work we investigate the relation between these S'-functions 

and the Lie symmetries of a rational 20DE. Based on this relation we can construct 
(^ ' a semi-algorithmic method to find the Lie symmetries of a 20DE even in the case 

^^ ■ where it presents no Lie point symmetries. 
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1 Introduction 

Lie group analysis is perhaps the most powerful tool to study differential equations 
(DEs) [2, 3, 4, 5, 6, 7, 8, 9, 10]. S. Lie [2] demonstrated that the majority of the 
technics of solving differential equations could be unified under the same theoretical 
background: the invariance of the DEs being solved under a continuous group of 
transformations (a Lie group). Since its appearance in the end of the 19th century, 
the Lie analysis of DEs has suffered a fantastic grown up especially in the last 
decades. One of the main reasons is, probably, the appearance of the computer in 
the second half of the 20th century, since it turns long and complicated symbolic 
calculations into the simple pushing of a button. Unfortunately, in real life, nothing 
is that simple. For the computer have the job done, we need to 'tell' it exactly what 
to do. In other words, it is necessary to furnish an algorithm, i.e., a finite sequence 
of determined steps. However, in Lies method, there is no algorithmic procedure 
to solve the determining equations for the infinitesimals (i.e., the coefficients of the 
symmetry generators)^. The things go worse when we are dealing with an ODE not 
possessing Lie point symmetries (presenting only dynamical symmetries). In this 
case we can not separate the determining equation in the derivatives and, instead 
of determining equations, we have only one determining equation: a 'very ugly' 
partial differential equation (PDE) that may leave us completely lost. In this last 
situation, we can not count even with a sistematic way to search for the symmetries; 
and, witout the symmetries, the Lie method can not be applied. 

To overcome the difficulties in the treatment of ODEs (or systems of ODEs) 
that do not possess Lie point symmetries several approaches have been developed: 
P. J. Olver introduced the concept of exponencial vector field (see [3], p. 185); 
B. Abraham-Shrauner, A. Guo, K.S. Govinder, P.G.L. Leach, F.M. Mahomed, 
A. A. Adam and others worked with the concept of hiden and non local symetries 
[13, 14, 15, 16, 17, 18]. C. Muriel and J.L. Romero have developed the concept 
of X-symmetry [19, 20] and E. Pucci and G. Saccomandi created the concept of 
telescopical symmetry [21]. Another great approach was brought by M.C. Nucci by 
making use of the Jacobi last multiplier [22]. 

Despite all these developments, there is still no fully algorithmic method to solve 
the determining equations. This mean that, in some part of the whole process, we 
can face a set of PDEs (the determining equations) for the infinitesimals we don't 
know how to deal with. 

In [1] we have developed an extension of the Prelle-Singer method [23] and, 
in that paper, we have proposed to use an unknown function (that we called S 
^) in order to make the 1-form^ (j)(x,y,y') dx — dy' proportional to an exact 1- 
form. In [27], we constructed a semi-algorithm to determine the S'-function for an 
20DE presenting an elementary^ first integral. In this paper we study the relation 
between the S'-functions and the Lie symmetries of a rational 20DE. Based on this 
relation we propose a semi-algorithm to calculate the Lie symmetries for rational 
20DEs. This procedure can succeed even in the case where there are only dynamical 



^Many computer packages implement some heuristics to find the infinitesimals and apply the Lie 
method. See, for example, [11, 12]. 

"^This idea was further pursued by [24, 25, 26]. 

'''This 1-form is associated with the rational 20DE y" — (p{x,y,y'), where is a rational function of 
{x,y,y'). 

^For a formal definition of elementary function, please see [28] . 



symmetries (i.e., no Lie point symmetries). 

In this work the main idea is to present a connection between a Darboux type 
method and the search for Lie symmetries of 20DEs. Here, we present one such 
algorithm that we expect is the firt of a series of developments in that area. 

The paper is organized as follows: In section 2, we present some basic definitions 
and stablish the relation between the S'-functions and the Lie symmetries of a 20DE. 
In section 3, we propose a method (a semi-algorithm) to calculate a Lie symmetry 
associated with an S'-function. In section 4, we present two examples to show the 
method in action. Finally, we present our conclusions and point out some directions 
to further our work. 

2 Lie Symmetries and /S-functions 

In this section, we will describe how the 5-function is defined and the relation 
between the Lie symmetries and the S'-functions of a 20DE. 

2.1 The definition of the 5- functions 

Let us consider the 20DE given by: 

y" = <t>{x,y,y'), (1) 

where is a function of {x,y,y'). 

A function /(x, y, y') defines a first integral (conserved quantity) of (1) if /(x, y, y') 
is constant over all solution curves of (1). So, in other words, the 1-form defined by 
uj = dl = Ix dx + ly dy + lyi dy' is null if the 1-forms (a, /3) defined by: 

a = (t)dx — dy' 

(3 = y'dx-dy (2) 

are null. This implies that 

UJ = dl = ra + s(3 (3) 

where r and s are functions of {x,y,y'). Thus 

dl = Ixdx + lydy + lyidy = r{(pdx — dy') + s{y'dx — dy) (4) 

implying that Ix = r (/) + s y' , ly = —s and ly' = —r. Therefore, if we determine r 
and s, we can find / via quadratures (see [29, 30]). 

Let us make some definitions that make it easier to show some connection be- 
tween the coeficients r and s and the symmetries of the 20DE (1). 

Definition 1: A function R{x,y,y') satisfying 

R{Adx + Bdy + Cdy') = d-f, (5) 

where A, B and C are functions of {x,y,y') and dj is an exact 1-form, is called 
an integrating factor for the 1-form {Adx + B dy + C dy') . 



Definition 2: Let S{x,y,y') be a function defined by 

S^-, (6) 

r 

where r and s are functions satisfying (3). We will call it a S -function associated 
with the 20 DE (1). 

Rewriting (4), we have 

dl = Ixdx + lydy + ly'dy = r ( (</> + - y')dx dy — dy' | . (7) 

We can see that r is an integrating factor for the 1-form ((0 + - y')dx — - dy — dy') . 
So, writing r = R and using the definition for the S'-function we can finally write 

dl = R {{(j) + S y')dx - S dy - dy') . (8) 

2.2 The connection of the S'-functions with the Lie sym- 
metries of a 20DE 

We will begin this section by stating a theorem (our first goal): 

Theorem 1: Let y" = (pix, y, y') be a 20DE presenting a first integral L{x, y, y') . 
If ri{x, y, y') is the infinitesimal of a Lie symmetry generator in the evolutionary 
form, then the function S defined by 

S^-^M, (9) 

Tl 

where D^ is the operator defined by 

Dx = dx + y' dy + (l){x, y, y') dy' , (10) 

is a S-function associated with the 20DE y" = (j){x, y, y') . 
To prove this result we will use the following lemma: 

Lemma 1: Let y" = (f){x,y,y') be a 20DE. If ri{x,y,y') is the infinitesimal of 
a Lie symmetry generator in the evolutionary form, then r] must obey the following 
PDE: 

Dl[ri]=DM<t>y'+^(t^y, (H) 

where D^ is the operator defined by (10). 

Proof of Lemma 1: If the hypothesis of the lemma is fulfilled, then 

X(2)[/-0(x,y,y')]=O, (12) 

where X^"^' = ^dx + rjdy + rj^^' dy' + rj^^' dy" is the second extension of the group 
generator X = S^dx + rjdy and rj^^' and r/'^^ are given by 

7?W ^ Dx[r]]-y'Dx[C], (13) 

^(2) = ^^[^(1)] _ y'axi^^ = ax [DxM - y'DxlC]] - y'Dxi^]. (14) 



Since X is a symmetry generator for the 20DE y" = (j){x, y, y') then the vector 
field XW defined by 

= C dx + ri dy + v'^^'' dy' - pdx - py' dy - pcpdyf = 

= {^-p)Ox + {v- py') dy + (??(^) -p4>) dy: = 

= (e - p) a, + (7? - y'p) dy + {Dx[n - y'p] - y' Dx[i - p]) dy,, (15) 

is also a symmetry generator. Choosing p = ^ we obtain a symmetry generator 
X in the evolutionary form: 

X'^^^ =r]dy + DMdy', (16) 

where r] = rj — y'S, ■ So, we have that X =ridy + D^ [^] dyi + D'^ [rj] dy// and we can 
write (12) as 

Dl[rj]-DM<t>y'-V(l>y = 0, (17) 

as we want to demonstrate. □ 

Now we can prove theorem 1: 

Proof of Theorem 1: If the hypothesis of the theorem is fulfilled, then we have 
that (see section 2.1, eq.(8)) 

dl = R{{(t>+ S y')dx - S dy - dy') . 

So, we can write Ix = R{4> + S y'), ly = —S R, lyi = —R. Using the compatibility 
conditions {Ixy — lyx = 0, Ixy' — ly'x = and lyyi — ly'y = 0), we get: 

Ry{cl) + Sy') + R{(l)y + Syy') + {SxR + SRx) = 0, (18) 

Ry> {4> + Sy') + R {<Py> + Sy:y' + S) + Rx = 0, (19) 

-{Ry:S + RSy:)+Ry = 0. (20) 

Eq.(19) plus eq.(20) times y' results 

Rx + y'Ry + (t>Ry'+R{4>y' + S)=0, (21) 

and eq.(18) minus eq.(20) times (j) results 

S {Rx + y'Ry + (f>Ry>) + R{Sx + y'Sy + 4>Sy') + R4>y = 0. (22) 

These equations can be written, respectively, as 

Dx[R] + R{4>y' + S) = 0, (23) 

SDx[R] + RDx[S] + R(py = 0. (24) 



where Dx is the operator defined in (10). Isolating D^ [R\ in eq.(23) and substi- 
tuting in eq.(24) we have that the 5-function must obey the following equation: 

Dx[S]=S'^+(Py,S-(l)y. (25) 



Now, let 17 be an infinitesimal defining a Lie symmetry in the evolutionary form. 



Substituting S = =^ in (25) we have 



.,^ .f^l^2 / n r^ix n2r^- 

D 



DM 



V 



DAv]\\^ f DM\ . _ DIM _ ^ DM 



_ . V \ — ] Yy ^ — — 'ry' _ Yy • 

1] J \ 7] J 7] 7] 



(26) 

By using lemma 1 we can verify that eq.(26) is an identity and the theorem is 
demonstrated. □ 

If the reader is familiar with the classical theory of ODEs he/she can note two 
interesting facts: 



Since the operator D^ represents (over the solutions of the 20DE) the total 

A) 



derivative with respect to x (i.e., -j-), the PDE for the S'-function is formally 



a Riccati ODE^. 

• If we apply the transformation u = —v'/v into the Riccati ODE u' = u'^ + 
g{x) u — h{x) we get the linear 20DE given by v" = g{x) v' + h{x) v. This 
20DE, is formally analogous to the linear second order PDE (linear 2PDE) 
for rj. 

In the next section we will use the relation (9) to construct a semi-algorithmic 
method to calculate the Lie symmetries of a rational 20DE. This method can be 
applied even to the case where there are only dynamical symmetries. 

3 A method to calculate the Lie symmetries 

The procedure can be divided in two main parts: first, we use the method 
developed in [27] to calculate the 5-functions associated with the 20DE; then we 
use the relation (9) to calculate the symmetries. 

3.1 Calculating the ^-functions 

Based on theorem 2 of [27] we have the following result''': 

Theorem 2: Let 

y =(t>[x,y,y) = —- -, (27 

N{x,y,y') 

be a rational 20DE (i.e., M and N are polynomial functions of {x,y,y')), and 
Let L{x, y, y') he an elementary first integral of it. Then, there is a S-function (S) 
associated with the 20DE (27) such that: 

(i) S is an algebraic function of {x,y,y'). 



®A Riccati ODE can be generally written as u' — f{x) u^ + g{x) u + h{x) 
"^For a proof, please see [27] 



(ii) The polynomial that defines the algebraic function S is an eigenpolyno- 
niial of the operator 

V = {N)D + (N^ 52 + {Ndy>M -M dy>N) S + N dyM - M dyN) dg, 
where D = N D^ (D^ is the operator defined by (10)). 

The claiming [ii) is simply a more formal way of saying that the S'-function must 
obey (25). Besides, it turns clearer what to do to obtain the S'-function: we 'only' 
have to calculate the eigenpolynomials (the Darboux polynomials) of the operator 
T). Each eigenpolynomial defines an algebraic function of (x, y, y') which is a S'- 
function associated with the 20DE (27). The emphasis on the word only is just to 
remember that calculating the eigenpolynomials of the operator T) may not be an 
easy task. But, although the procedure may be hard to apply, it is of an algorithmic 
nature. 

In this paper we want to deal with 20DEs integrable by quadratures (i.e., pre- 
senting Liouvillian^ first integrals). Since we do not know the general form of the 
S-funtions for this case, we will restrict ourselves to the case where the 20DE 
presents rational S-funtions. We will show that even with this restriction, we can 
find symmetries in a lot of interesting cases. 

3.2 Finding the symmetries 

Once we have found the rational S- functions associated with the 20DE (27), 
we will use the relation Dx^]/ri = — S to obtain a Lie symmetry. As we have 
mentioned in the last section, we will restrict ourselves to the case where the 20DE 
(27) presents rational S-funtions. For this case we were able to prove the following 
result (this paper's second goal): 

Theorem 3: Let y" = (^(x, y, y') = M/N (M and N polynomials of (x, y, y') ) 
be a rational 20DE presenting two independent Liouvillian first integrals Ii(x,y,y') 
and l2{x,y,y') . If this 20DE has two rational S-functions Si and S2 associated 
with it such that dli = Ri{{(t>+ Si y')dx — Si dy — dy') and dl2 = i?2((0+ S2 y')dx — 
S2 dy — dy'), then the following statements hold: 

(i) There exists two independent Lie symmetries in the evolutionary form 
(16) such that the infinitesimals rji and 7/2 are Darboux functions of {x,y,y') , 
i.e., the infinitesimals r/i and 1^2 have the form 

m = e^i UpuT, (28) 

i 

V2 = e^ Up^T- (29) 

j 

where Ai, Bi, A2, B2, the pn and the p2j are polynomial functions of {x,y,y') and 
the c\i and C2j are constants. 

(ii) The pn and the p2j are irreductible eigen-polynomials (Darboux poly- 
nomials) of the operator D (defined by N Dx ) or are irreductible factors of the 
denominators of Si, S2 , respectively. 



*For a formal definition of Liouvillian function, please see [28]. 



(in) The irreductible factors of the polynomials Bi and B2 are irreductible 
eig en-polynomials (Darhoux polynomials) of the operator D or are irreductible fac- 
tors of the denominators of Si, S2 , respectively. 

In order to prove this theorem we will need some results: 

Proposition 1: Let y" = (j){x,y,y') be a 20DE presenting a first integral 
I{x,y,y') . If S{x,y,y') is a S -function associated with it, then the lODE defined 
by 

= -S{ai,u,v) (30) 

du 

has I{ai,u,v) = C (where C is a constant) as its general solution. 

Proof of Proposition 1: I{ai,u,v) = C is a general solution of the lODE 
dv/du = —S{ai, n, v) if and only if 

Du[I] = 0, (31) 

where D,^ = du — S d^. The condition (31) means that dul — Sd^I = 0, i.e., 
S = lu/Iv ■ Consider that the hypothesis of proposition 1 are satisfied. To prove 
proposition 1 we only have to prove that S = ly/Iy' ■ According to (8) we can 
write dl = R{{(j) + S y')dx — S dy — dy'), leading to ly = —R S and ly/ = —R . So, 

S = ly/Iy,. U 

Using similar reasoning we can enunciate the following^: 

Proposition 2: Let y" = (j){x, y, y') be a 20DE presenting a first integral 
I{x,y,y') . If S{x,y,y') is a S -function associated with it, then the lODE defined 

by 

dv 

-— = (j){t,a2,v)-\-vS(t,a2,v) (32) 

has I{t,a2,v) = C (where C is a constant) as its general solution. 

Proposition 3: Let y" = (j){x, y, y') be a 20DE presenting a first integral 
I{x,y,y') . If S{x,y,y') is a S -function associated with it, then the lODE defined 
by 

du ^ 4>(t,u,a3) +a3S{t,u,a3) 

dt S{t,u,a3) 

has I{t,u,a3) = C (where C is a constant) as its general solution. 

Definition 3: The three lODEs (30), (32) and (33) will be called auxiliary 
lODEs of the 20DE (27) associated with the S-function S{x, y, y') . 

In what follows, we will refer to (30) as lODE Al, to (32) as lODE A2 and to 
(33) as lODE A3. This concept (of auxiliary lODE) will be useful to enunciate 
the following result. Note that the general solution of the auxiliary lODEs (Al, 
A2 and A3) are defined by the function / that is a first integral of the 20DE. 
The following result establishes another important link between the 20DE and the 
auxiliary lODEs. 



'The proofs of propositions 2 and 3 are analogous to that of proposition 1. 



Proposition 4: Let y" = (j){x,y,y') be a 20DE presenting a Lie symmetry that 
can be written (in the evolutionary form) as 

X^'^ =ri{x,y,y')dy + D, [rj(x, y, y')] dy, . (34) 

Let C^ = Dx\ff\ , then the operator Y = r](ai,u,v) du + Cio-iiUjv) d^ is a generator 
for a Lie point symmetry of an auxiliary lODE Al (30). 

Proof of Proposition 4: The condition for a vector field Y = ri{ai,u,v) d^ + 
C{ai,u, v) dv to be a symmetry vector of the lODE (30) is that its commutator with 
the operator Du = du — S dy is proportional to D^- So we have to prove that [Y, Z)„] 
is proportional to D^- We have that: 

[Y, Du] = -ivSu + CSv) d, - i^u - Sv.) du - (Cu - SU d, = 

= -Du[v]du-{riSu + CS, + Cu-SUdu. (35) 

Since S = —Dx[ri]/ri = —Q/rj we can write 

\ T T J 

= -Du[ri]du-{-SDu[v])d,, (36) 

implying that [y,D„] = -!)„[?/] D„. □ 

Proposition 5: Let y" = (f)(x,y,y') be a 20DE presenting two independent 
first integrals Li{x,y,y') and l2{x,y,y') . If Si = Liy/Li^ and 5*2 = I2y/l2z are 
two S'-functions associated with it such that dLi = Ri{{4> + "S*! y')dx — Sidy — dy') 
and dL2 = R2{{4> + S2y')dx — 5*2 dy — dy'), then there are infinitesimals for Lie 
symmetries in the evolutionary form, rj;^ and % , (such that —I^x^i\Mi = Si and 
-Dx\ri2]M2 = S2 ) given by 

"^^ = {Sl-S2)R2' ^^^^ 

'^'= {S2-'si)Ri - ('') 



Proof of Proposition 5: Consider that the hypothesis of the proposition are 
fulfilled. Let the operators Xi =riidy + Dx[rii] dy' and X2 = 7/2 5^ + -Da;[%] dyi be 
Lie symmetries of the 20DE y" = (pix, y, y') in the evolutionary form. Without 
loss of generality we can choose r]i and 7/2 such that ^i[/i] = -^2 [-^2] = and 
Xi[l2] = X2[Li] = 1 . Using Xi[Li] = and Xi[l2] = 1 we have 

fjl Ily + Dx[ril\ Llyl = 0, 

m hy + D^iVi] hy' = 1, 



implying that 



-Ih 



^1 = -T-T -J—T = -1 = —. TT^ ■ (39) 

hy hy' - hy' hy ^ l2y, - l2y h 



Ly, ■' ~ a \^^y' '2y, 



Since l2y' = —R2 , hy/hy' = Si and I2y/l2y' = 5*2 we can write rji = l/{{Si — 
S2) R2) and we can deduce the format (38) for r/2 in an anafogous way. □ 

Theorem 4: // a rational lODE written in the form y' = M{x,y)/N{x,y) 
(M and N are polynomials with no common factor) has a general solution of the 
form I{x, y) = C where I is a Liouvillian function of its arguments and C is an 
arbitrary constant, then the 1-form M dx — N dy presents an integrating factor R 
of the form 

i^ = e'-«^^^'^ (40) 

i 

where tq is a rational function of {x,y), the pi and the factors of the denominator 
of ro are irreductible Darboux polynomials of the operator N dx + M dy and the Ci 
are constants. 

For a proof of theorem 4 see [31, 32, 33, 34] or [35, 36]. 

Now we can prove theorem 3. The key point in the fohowing demonstration 
is the fact that the general solutions of the auxiliary lODEs associated with S 
are defined by the same function / that is a first integral of the rational 20DE 

y" = <P{x,y,y')=M/N. 

Proof of Theorem 3: First of all let's stablish some notation: we will write 
^i = Pi/Qi and 5*2 = P2/Q2 , where Pi , Qi , P2 and Q2 are polynomials of 
(x, y, y') and the pairs Pi , Qi and P2 , Q2 do not have any common factors. 

Consider now that the hypothesis of the theorem are fulfilled. From proposition 
1, we have that /i(ai,n, f) = C is a general Liouvillian solution of the rational 
auxiliary lODE Al associated with Si 

dv „, . Pi(ai,u,w) 

— = -Si{ai,u,v) = -—- -. (41) 

du Qi[ai,u,v) 

Then (from theorem 4) the 1-form Pi(ai, n, v) du + Qi{ai,u, v) dv presents an inte- 
grating factor of the form (40). Let Pj^i .5i](ai,u, f) be such an integrating factor. 
Since 



(42) 



dh{x,y,y') 


= Pi ((0 + 5i y')dx -Sidy- dy') 






= -^i{MQi+PiNy')dx-PiNdy- 


-QiNdy'), 1 


we can see that 






dh , „ 


Ri{x,y,y')Pi{x,y,y') ^ dh , 
- Qiix,y,y') '^^ Oy'^^^y^y^- 


= -Ri{x,y,y') 


implying that 








Ri{ai,u,v)Pi{ai,u,v) . dh , . 

and (ai,n,v = 
Ui[ai,u,v) ov 


= —Ri{ai,u,v). 



Therefore, —Ri{ai,u,v)/Qi{ai,u,v) is also an integrating factor for the 1-form 
Pi(ai, u, v) du + Qi{ai,u, v) dv . Then, we can write ^ = J-'i{h) R[Ai,si] (where 



10 



T\ is a function of /i ) or n^ = fei -^[^i,^!] (where k\ is a constant) leading to 

Rx=Tx{h)R\A\,s^\Q\ (or i?i = fei i?[Ai,si] <3i )• So, i?[Ai^5i](x,y,y') Qi(a^>y>y') 
is an integrating factor for the 1-form (</> + S\ y')dx — Sidy — dy' . 

Besides, from proposition 2 we have that Ii{t,a2,v) = C is a general Liouvillian 
solution of the rational auxiliary lODE A2 associated with Si 

dv \ , au \ Qi{t,a2,v)M{t,a2,v) + vPi{t,a2,v)N{t,a2,v) 

— = (f)[t,a2,v) + vbi{t,a2,v) — 



dt Qi{t,a2,v) N{t,a2,v) 

(43) 
Then (from theorem 4) the 1-form —(QiM + v Pi N) dt + QiNdv presents an 
integrating factor of the form (40). Let R[yi2,si](t,a2,v) be such an integrating 
factor. From (42) and using a completely analogous reasoning of what we have 
shown above, we can conclude that Ri{t,a2,v)/ {Qi{t,a2,v) N{t,a2,v)) is also an 
integrating factor for the 1-form —{Qi M+v Pi N) dt+Qi N dv . Then, we can write 
^Iv = -^2 (A) R[A2,si] (where T2 is a function of h ) or -^^ = k2 R[A2,si] (where 
^2 is a constant) leading to Ri = J^2ih) R'[A2,si] Qi ^ (o^ -^1 = ^2 R[A2,si] Qi -^ )• 
So, Riji2,si]{^iyiy')Qii^iyiy') ^i^iUiy') is an integrating factor for the 1-form 
{(j) + Si y')dx - Sidy - dy' . 

Now let's see: since R[Ai,si]{x,y,y')Qi{x,y,y') and R[A2,si]{x,y,y') Qi{x,y,y') 
N{x, y, y') are integrating factors for the 1-form {(f)+ Si y')dx — Sidy — dy' , then we 
can write R[ai,si\ {x, y, y')Qi{x, y, y') = F{Ii) R[A2,si\ {x, y, y') Qi{x, y, y')N{x, y, y') 
( J" is a function of h ) or R[ai,si\ {x,y,y')Qi{x,y,y') = ki R[A2,si] {x,y,y') Qi{x,y, 
y') N{x,y,y') {k is a constant). 

• First possibility: R[Ai,si]{x,y,y') = k R[A2,si]{x,y,y') N{x,y,y') . Let's re- 
member that the polynomials forming R[Ai,si]ix,y,y') are polynomials of 
{y, y') and the polynomials forming R\a2,si] {x-, y-, 2/') are polynomials of (x, y') . 
Therefore, we can conclude that R[Ai,si]{x,y,y') and R[A2,si]{x,y,y') are of 
the form exp[ro] Hi Pi'^* 1 where the pi are polynomials of {x,y,y') and vq is 
a rational function of {x,y,y') . 

• Second possibility: R[jii^si]{x,y,y') = J^ih) R[A2,si](.x,y,y') N{x,y,y') . For 
the case where Ii is an elementary function of (x, y, y') , from the results 
presented in [29] we can conclude that R[Ai,si]{x-,y-,y') and R[A2,si]{x-,y-,y') 
are of the form Hi Pi^^ and use the same reasoning above to show that the pi 
are polynomials of (x, y, y') . 

Now, let's suppose that Ii is a non elementary Liouvillian function of (x, y, y') . 
Since R[jii^si]{x,y,y') and R[A2,si]ix,y,y') are elementary functions, the only 
possibility for the relation R[ai,si] {x, y, y') = ^{h) R[A2,si] {x, y, y') N{x, y, y') 
to be true is that J'{Ii) is the constant function, i.e., J^{Ii) = k. So, 
we are again in the first possibility and we have that ^mi,51] (2^1 ?/;?/') and 
R[A2,si]ix,y,y') are of the form exp[ro] Hi Pi"^^ ) where the pi are polynomi- 
als of (x, y, y') and tq is a rational function of (x, y, y') . 

Now, consider the auxiliary rational lODEs Al and A2 associated with S2 '■ 

dv P2{ai,u,v) 

— = -S2{ai,u,v) = ——- -, (44) 

du Q2[ai,u,vj 
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dv \ , ou \ Q2it,a2,v)M{t,a2,v) + vP2{t,a2,v)N{t,a2,v) 

— = 4>{t,a2,v) + vS2{t,a2,v) = — — — — r . 

at Q2{t,a2,vjN{t,a2,v) 

(45) 
From theorem 4 (and using propositions 2 and 3 again) , we can infer that the 1-forms 
P2{ai,u,v) du+Q2{ai,u,v) dv and -{Q2{t,a2,v) M{t,a2,v)+v P2{t,a2,v) N{t,a2,v 
)) dt + Q2(t,a2,v) N(t,a2,v) dv present integrating factors of the form (40). Let 
R[Ai,S2]iO'i,u,v) and R[A2,S2]{t,CL2,v) be, respectively, these integrating factors. 
Using the same reasoning above we can conclude that (analogously as above): 

• ~R\Ai,S2]{^iyiy')Q2{x,y,y') is an integrating factor for the 1-form ((/> + 
S2 y')dx — S2dy — dy' . 

• R\^A2,S2]{x-,y-,y')Q2{x-,y-,y') N{x,y,y') is an integrating factor for the 1-form 
{4> + S2 y')dx - S2dy - dy' . 

• R[Ai,S2\{x,y,y') and R[A2,S2\{x,y,y') are of the form exp[ro] Hi Pi'^S where 
the Pi are polynomials of (x, y, y') and rg is a rational function of (x, y, y') . 

Therefore, from these results and the analogous above, we can afirm that there 
exists two independent Liouvillian first integrals Ii and I2 such that 

dli{x, y, y') = Ri {{(t> + Si y')dx -Sidy- dy') , (46) 

d72{x, y, y') = R2 ((</> + S2 y')dx - S2dy - dy') , (47) 



/^ 



where Ri = —R[ai,si] Qi and R2 = —R[ai,S2] Q2 are of the form e '^ Hi P- 
{A,B,pi polynomials in {x,y,y') and Cj constants). 

We can now prove (i). From (46), (47) and proposition 5 we can conclude (since 
5*1 and S2 are rational functions of {x,y,y')) that there are infinitesimals rjj^ and 
rJ2 of the form (28) and (29). 

In order to prove (ii) and (iii) let's first note that 

D^n ^ Q^m = _^p,, (48) 

vi Qi vi 

-BM = ^ =, Q2£M = _NP2, (49) 

V2 Q2 V2 

i.e., Qi -'^' and Q2 - are polynomials in {x,y,y') . So, we have the situation: 

the infinitesimals fji and 7^2 are Darboux functions such that -''-' = poll and 

-^' = P0I2 , where Di = Qi D and D2 = Q2 D ■ Finally, from the main result 
presented in [33], we can directly conclude (ii) and (iii) . □ 



The procedure to calculate the symmetries in this case is completely analogous to 
that for calculate the integrating factors of rational lODEs with Liouvillian solution 
(see [34]): briefly, we calculate the Darboux polynomials (up to a certain degree) of 
the operator D . Then we substitute the candidates for rj in the equation 

DM 
Q^ = -NP, (50) 

V 

and we solve the resulting linear equations for the constants Cj . 
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4 Examples 



In this section we will present two examples where we can use the semi-algorithm 
presented above to calculate the Lie symmetries. In the first example we will show a 
20DE with one point symmetry and one Darboux dynamical symmetry to ilustrate 
theorem 3. This symmetry can be found by our method - a process devoided of 
any guess. We will use this example to show our procedure in action, i.e., we will 
provide comments about the calculations to clearer the steps of the method. In the 
second example we use a 20DE not presenting Lie point symmetries but with two 
rational dynamical symmetries. The good point to be noted is that, for our kind of 
process, the existence or not of point symmetries makes no difference. 

4.1 First example 

Consider the 20DE given by 

„ 2y -'izy + z^y- ZX + z^x 
y = -( ^ • (51) 

To use the procedure explained in the section 3 we have, first, to obtain the operators 
D and T> . They are given by 

D = {xy -y^)— + (xyz-y^z)— + {3yz-yz^ - xz^ -2y + xz) — , (52) 

d d 

V = {y'^ + y'^x^ - 2y^x)—- + {-2zy^x + z/ + zy^x^)—+ 

ox oy 

{yzx^ - yx^z'^ + 2^ + 2 y'^zx + y^z^ -3zy^ -2 xy'^)—+ 

oz 

({y'^x^ + / - 2 y^x) s^ + (-2 zyx"^ + 2zy^ + 2 xy^ + yx'^ - 3 y^) s 

-z^x^ + zx^ + z^y^ - 2 zyx + 2 y^ - 3 zy^ + 2 z^yx\ —- . (53) 

^ OS 

The second step is to calculate the Darboux polynomials of the operator V . We 
can find two Darboux polynomials: 

Psi = {z-l){y + x) + sy{x -y), (54) 

Ps2 = x{-2y + zy + zx){z - 1) + sy{x - y){zx - y). (55) 

The S'-functions are the solutions of Darboux polynomial of T> = 0. From (54) and 
(55) we can obtain two rational S'-functions given by 

{z-l){y + x) 

51 = ^ — , (56 

y{x-y) 

g2 = x{-2y + zy + zx){z-l) 
{x - y) (zx -y)y 

The next step is to calculate the Darboux polynomials of the D operator and the 
corresponding cofactors to build the infinitesimals. They are 

Pi = y =^ qi = {x - y) z, (58) 
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P2=x-y =^ q2 = -(z-l) y, (59) 

P3 = z-1 =^ q3 = 2y - zy - zx. (60) 

Using these Darboux polynomials and looking at (50), we can use the semi-algorithm 
described in [34, 37] to calculate the infinitesimals: 

{x-yfe(-^^y 
Vi = , (61 

y 

{zx -y){x- yf 
V2 = 7 7\ • (62 

{z-l)y 



We can use the relations ^i[/i] = , X2[Ii] = 1 , D[Ii] = and solve them to 
Iix , hy and /i^ . We can also use the relations Xi[/2] = 1 , ^2(^2] = , D[l2] = 
and solve them to l2x , l2y and l2z ■ From the derivatives we can integrate and 
obtain /i and I2 '■ 

h = ^^. (63) 

[y-x) 

revi^'-i) [-2y + zy + zx) y 
l2 = j ^— p dx. (64) 



4.2 second example 

It can be proved (see [19]) that the 20DE 

„ x^ + iy^ + 2y^ 

y = 4? — (6^) 

has no Lie point symmetries. By doing the procedure explained in the section 3 we 

have: 

The operator D and D are given by 

D = 4y^^ +Ayh^ + (-x' -4y^ -2y^) ^, (66) 

ox oy ^ ' oz 

V = 16/— + 16y6^— + 4y2 (_4y5 -2y^ - yx^) — + 
ox oy ^ ' oz 

+4y2(-3x2+4yV-2y2+4/)— . (67) 

We can find two Darboux polynomials of the operator D : 

Psi = x + yz + Sy'^, (68) 

p^2 = 4y^x + 8y^z + 4:y'^x + 4y^z + x^ + 2 x^yz + 2 i {x^y^ + 4 /)+ 

+4Sy^[-y + zx + 2yz'^ + 2iy'^z). (69) 
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From (68) and (69) we can obtain two rational S'-functions given by 

SI = -^. (70) 



A {—y + zx + 2yz'^ + 2iy'^z)y^ 
The Darboux polynomials of the D operator and the corresponding cofactors are 



Pi=y =^ qi=Ay^z, (72) 

P2 = x + 2yz + 2iy'^ =^ q2 = -2y{x - 2yz - 2iy'^), (73) 

ps = x + 2yz-2iy^ =^ qs = -2y{x - 2yz + 2iy^). (74) 

Using these Darboux polynomials, we can use (50) to calculate the infinitesimals 

_ ^ V 

^^ {x + 2yz-2iy^){x + 2yz + 2iy'^)' ^ ^ 

_ ^ _ -y + zx + 2yz'^ + 2iy'^z 

"^2- A{x + 2yz + 2iy^) ' ^ ^ 

We can use the relations Xi[Ii] = , X2[h] = 1 , D[Ii] = 0, Xi[/2] = f , 
X2[l2] = , D[l2] = to obtain h and h : 

/,=4.-2nnfi±|i^i^y (77) 

yx + 2yz + 2ty^ J 

h = -y"^ + -ix + ln{y) - -Inix + 2yz + 2iy^) - -^ + -z\ (78) 

In [19] Muriel and Romero were capable of dealing with this lODE (65) by 
developing the concept of C°°(M(^')-symmetries - vector fields (also called A- 
symmetries) that are neither Lie symmetries nor Lie-Backlund symmetries. The 
determining equations for the components of these vector fields depend on an arbi- 
trary function A , which can be chosen in order to simplify the process of solution of 
the determining equations. For this lODE they obtained A = x/u^ and v = x-^ 
for the A-symmetry. 

In our process (shown above) the main cost (computational cost) is to calculate 
the Darboux polynomials of the operator T) (followed by calculating the Darboux 
polynomials of the operator D). The good news are that the process is entirely 
computational, i.e., devoid of any guess. 

5 Conclusion 

Although the Lie symmetry method is the most powerful method for solving, reduc- 
ing and studying dynamical systems since its appearance (end of XIX* century), 
there are (still today) a bunch of open questions. The main 'gap' is the absence of 
an algorithm to calculate the symmetry vector fields. The things go worse when 
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the dynamical system under study does not present Lie point symmetries (in which 
case we can not even count with a systematic method to deal with the problem). 

In this paper we have proposed a semi-algorithm to find Lie symmetries for a 
rational 20DE. We began by showing a deep connection between the S'-functions 
(see section 2.1 and [1, 27]) of a 20DE and its symmetries in the evolutionary form. 
Then, restricting ourselves to the case where the S'-functions are rational and the 
first integrals are Liouvillian, we could develop a semi-algorithm to calculate the 
symmetries. 

The great advantage of our method is that it converts the search for symmetries 
into (essentially) searching the Darboux polynomials of the polynomial differential 
linear operators V and D (in other words, into solving second degree algebraic 
equations) - a semi-algorithmic procedure. The disadvantage is that, as the degree 
of the Darboux polynomials grows, the computational coast increases 'a lot' and, 
sometimes, turns the entire process inviable. 

Our method is not general since it is limited to the cases where the S'-functions 
are rational (see section 3). However, from a practical point of view, this restriction 
does not seems to be a great one: all the non-linear rational 20DEs (that we have 
analysed) with two Liouvillian first integrals, presented rational S- functions. 

In regard to future work, there are some main directions (open questions): 

• What is the general form of the S-functions for a 20DE presenting two Liou- 
villian first integrals? (and only one Liouvillian first integral?) 

• What about NODEs? 

• What happens if we consider elementary functions present in the ODE? 

• What about systems of ODEs? 

• Etc... 
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